Magnetic phase diagram of the Kondo lattice model with quantum localized spins 
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The magnetic phase diagram of the ferromagnetic Kondo lattice model is determined at T=0 in 
ID, 2D, and 3D for various magnitudes of the quantum mechanical localized spins ranging from S — 
i to classical spins (5* ^ cxd). We consider the ferromagnetic phase, the paramagnetic phase, and 



the ferromagnetic/antiferromagnetic phase separated regime. There is no significant influence of the 
spin quantum number on the phase boundaries except for the case S = | , where the model exhibits 
an instability of the ferromagnetic phase with respect to spin disorder. Our results give support, at 
least as far as the low temperature magnetic properties are concerned, to the classical treatment of 
the S = 3/2— spins in the intensively investigated manganites, for which the ferromagnetic Kondo- 
lattice model is generally employed to account for magnetism. 
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I. INTRODUCTION 

The ferromagnetic Kondo lattice model (FKLM), also 
known as s-d model or double exchange model, has drawn 
a lot of attention over the past years in the field of mag- 
netism and electronic correlations. The model consists 
of Bloch electrons coupled to localized spins sitting on 
the sites of a crystal lattice. For the case of strong 
(Hund) coupling and an energetically favored parallel ori- 
entation of a localized moment and an electron, Zener 
proposed the double exchange mechanism to explain fer- 
romagnetism (FM) in manganitesi The gain in kinetic 
energy of the conduction electrons favors a parallel con- 
figuration of the localized spins. In the framework of 
a two-site model, Anderson and Hasegawa showed that 
the hopping amplitude of the electrons between sites i 
and j is proportional to cos{9ij /2), where 6ij is the angle 
between the localized spins. ^ 

A major field of application for the FKLM is linked 
to the phenomenon of colossal magnetoresistancc^ in the 
manganese compounds already aimed at by Zencr. Here, 
the 5 Mn d-shells are split by the crystal field into three 
degenerate t2g-orbitals forming localized spins of 5 = |. 
They interact via Hund's rule with itinerant electrons 
stemming from the remaining two degenerate e^-orbitals. 
With a Hund exchange interaction mostly estimated to 
be several times the hopping amplitude'* ''''^ the mangan- 
ites belong to the rather strongly coupled materials. Al- 
though there are other important aspects to take into ac- 
count when modelling the rich physics of the manganites, 
like eletron-phonon-coupling and the second conduction 
band, the FKLM already in its simplest single band ver- 
sion is crucial for understanding at least the magnetic 
properties of this class of substancesf^ 

A hot topic where the FKLM is used as a basic model 
are the diluted magnetic semiconductors (DMS) with 
promising technical applications for microelectronics.^'® 
These materials consist of a (often HI/V as, e.g., GaAs) 
semiconducting host and substituted transition metal im- 
purities (e.g. Mn) occupying cation sites, the latters ex- 
hibiting ferromagnetism due to a coupling of the localized 
cation spins mediated by a spin exchange interaction with 
valence and impurity band holes. In the case of DMS, 



this interaction is considered intermediated 

There have been various types of approaches to solv- 
ing the many-body problem of the FKLM in order to 
get a phase diagram. On the theoretical side, several 
treatments are based on Dynamical Mean Field Theory 
(DMFT),i°i"d'^i" In Ref. [lla Schwinger-boson method 
is used and applied to 2D and 3D. Bosonization in ID 
is extensively discussed in Ref. |3 More recently an 
analytical continuum field theory in 2D was employed.^® 
Valuable information to compare theoretical results with 
can be gathered from Monte Carlo simulations >i2iiiii 
The main feature of the magnetic phase diagram is the 
same in all these works: with increasing coupling strength 
ferromagnetism (FM) prevails for all charge carrier den- 
sities except for a more or less small region around half- 
filling where antiferromagnetism (AFM) or FM/AFM 
phase separation exists. 

Most approaches rely on the assumption that the local 
spins can be treated classically. This simplification has 
been justified by checking the classical spin against quan- 
tum spin approaches, both giving similar resultsi^^ii^ 
More recently, a phase diagram in ID was obtained 
by means of the Density-Matrix Renormalization Group 
(DMRG) yielding numerically exact results for a quan- 
tum spin S — 1/2.^® However, we do not know of any 
systematic, quantitative analysis of the influence of the 
spin magnitude on the magnetic properties of the FKLM. 

In this work we present phase diagrams of the FKLM 
at T=0 for several spin quantum numbers in ID, 2D, and 
3D. We use an equation of motion approach exploiting 
exactly solvable limiting cases and exact relations among 
Green functions and among their spectral moments. By 
evaluating the (free) energy at T=0 we can distinguish 
three different phases: paramagnetic (PM), ferromag- 
netic, and ferromagnetic-antiferromagnetic phase sepa- 
rated (PS). The latter two are the typical phases in the 
strong coupling region of the phase diagram-^i2ii24i4iiSii& 
which is the regime relevant for the manganites and on 
which we want to focus in our work. It should be men- 
tioned that in the weak to intermediate coupling regime 
several other phases like canted, spiral, or island have 
been foundi^J^ 

One of the main results will be that there are no ma- 
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jor differences in the FM-PS phase boundaries for spin 
quantum numbers S > 1/2. The case S = 1/2 is spe- 
cial: here we obtain an instability of the ferromagnetic 
against the paramagnetic phase. To understand this be- 
havior we discuss the spectral weight distribution of the 
excitation spectrum and its modifications by a variation 
of the quantum character of the localized spins. 

The paper is organized as follows. After presenting our 
theoretical approach in section II we discuss the phase 
diagrams of the FKLM in ID, 2D, and 3D for different 
quantum numbers of the localized spins in section III. 
A summary and an outlook on possible improvements of 
our theory are given in section IV. 



II. THEORY 



The Hamiltonian of the FKLM reads 



H 



,t 



(1) 



The first term describes Bloch electrons of spin a with 
a nearest neighbor hopping integral t. c^^J annihilates 
(creates) an electron of spin a at lattice site i. The lat- 
tice is chosen to be simple cubic in our case, s is the 
electron spin and S the localized spin operator, and both 
are coupled through a Hund exchange J. Using standard 
second quantization notation the interaction term can be 
rewritten 



Hi, 



S. 



(2) 



with = ±1, ' arc the z-component, raising and 
lowering operators for a localized spin at site «, and ni„ = 
'^ia^iu is the occupation number operator at site i. 

The many-body problem of the above Hamiltonian is 
solved with the knowledge of the one-electron Green func- 
tion G\^a{E), or, equivalently, the electronic self-energy 
^^.a{E)■. 



k 

fl 



E - e(k) - Ska(£;) ' 



(3) 



with the Bloch dispersion e(k). One can then calculate 
the internal energy U of the FKLM, being equivalent to 
the free energy at T=0, for the ferromagnetic and the 
paramagnetic state. There is a simple relation between 
the energy U of the FKLM and the imaginary part of the 
corresponding one-particle Green function, 

1 /•+OC 

^ = = TJTiY. '^^f- (E)ESu. (E) , (4) 

where Siia{E) = —■^QGuciE) is the local spectral den- 
sity and /_ denotes the Fermi functioniiS 



Note that the existence of the ferromagnetic state is 
supposed and not the result of a self-consistent calcula- 
tion, i.e. the magnetization is a parameter in our scheme. 

The method wc chose to solve the Hamiltonian ^ for 
the Green function (PJ is a moment conserving decou- 
pling approach (MCDA), which does not require the lo- 
calized spins to be classical. This theory has been applied 
before in model studies^?, and to real substances^^''^^. For 
a detailed account of the decoupling procedure we refer 
the reader to Ref^. Here we summarize the main points 
of the method and want to emphasize features which are 
important for the following discussion of our results. 

The starting point is the equation of motion for the 
Green function ||2Jl. The generated higher Green func- 
tions read 



IikJcr{E) 
Fikja{E) 



as- 



'Ck-a\c]^))E ■ 



(5) 
(6) 



Iik,ja{E) is a Ising-like GF solely comprising the z- 
components of the spins, whereas FikjaiE) describes 
spin flip processes which are neglected when using classi- 
cal localized spins. After writing down the equations of 
motion of likjaiE) and Fikja-iE) the decoupling is per- 
formed. Of special importance for correlation effects is 
the treatment of the local higher Green functions, namely 



F^P.„{E) 



2) 
r(4) 



(E) 
(E) 



F^E) 



{{srs!c,-y,cl))E , (7) 

{{srs^c,.;cl))E, (8) 

{{Srn^.c,-y,cl))E , (9) 

((5fn,_,c,,;4))^. (10) 



These functions are expressed in terms of the lower order 
GF ©, lO, and © with coefficients fitted by the first 
two spectral moments of these GFs, respectively, repre- 
senting a non-perturbative approximation for the whole 
temperature rangc:^- The choice of the "correct" lower 
order GF is guided by some non-trivial limiting cases 
which we summarize next. 

For an assumed complete ferromagnetic polarization 
{{S^) = S) of the FKLM one obtains from the spectral 
representation of the Green functions: 



Fl-lm [[S-^]+^z.] FuME) ,(11) 



2j + 2^ 

fI^JE) ^'^=' SG,,AE) - z^IuME) ■ (12) 



For S = ■^j a case which is of particular importance to 
our investigation due to its maximum amount of quantum 
fluctuations, the following relations hold at any temper- 
ature (i.e. at any {Sz))- 



F^^iE) 'J \z„Fu,.{E) , 



(13) 



^G^,„{E)-z,Iu,ja{E) . (14) 
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Furthermore in the case of a fully occupied conduction 
band: 



n=2 



(15) 
(16) 



The above exact relations are used to motivate the 
following ansatz for the higher local GF: 

fS^(E) = ai.G.,.(i?)+/3i.F,,,,^(i?), (17) 

F^]AE) = a^.G^.^E) + p2.luME) , (18) 

F^]AE) - a3.G.,.(i?)+/33.F,,,,.(i?), (19) 

fI^^{E) - a^^G^j^E) + Pi„Iu,AE) ■ (20) 

The temperature dependent interpolation coefficients 
otim Pier (* = Ij ■•j4) depend on various correlation func- 
tions and are listed in Appendix A. It is easily verified 
that the approximations p7|) - (|20|l fulfill the exact limit- 
ing cases (I11II - H16II . In addition our approach reproduces 
the limit of the ferromagnetically saturated semiconduc- 
tor (T=0, band occupation n=0)i^ 

The resulting self-energy S(t(-B) is local and depends 
on various expectation values of pure fermionic, mixed 
fermionic-spin, and pure localized spin character: 



(S^), {{S^r), ((5^)3), {S+S-}) , 



(21) 



The site indices have been dropped due to translational 
invariance. Whereas the first two types can be calculated 
within the MCDA the localized spin correlation functions 
are known for ferromagnetic saturation and in the param- 
agnetic state (see Appendix B). The many-body problem 
represented by ^ can thus be solved approximately for 
the Green function in a self-consistent manner,— We 
emphasize that the quantum mechanical character of the 
localized spins is fully retained in our approach. Fur- 
thermore there is no restriction to the parameter range 
within which our method can be applied. 

In order to determine the phase boundary between the 
ferromagnetic and the ferromagnetic-antiferromagnetic 
phase separated region we first have to solve the Hamil- 
tonian for an antiferromagnetic configuration. Using 
the standard sublattice decomposition for bipartite lat- 
tices and neglecting the off-diagonal elements of the self- 
energy matrijiS^ one obtains the following Green function 
for sublattice A: 



GtiE) 



E-e'{k)~j:tiE) 



|t(k)P 



-(22) 



with the diagonal elements of the self-energy matrix 
Sj^^ = ^k-cr: the sublattice dispersion e'(k) and the 
inter-sublattice dispersion i(k). The approximate solu- 
tion for the self-energy presented above for the transla- 
tionally invariant case can be obtained analogously for 



the antiferromagnetic case. The energy of the antifer- 
romagnetic phase can be evaluated using by simply 
replacing Sua {E) by {E) . Averaging over the sublat- 
tices is not necessary due to symmetry reasons, i.e. the 
summation over the sublattices is absorbed into the spin 
summation. 

In this work we restrict our considerations concerning 
AFM to G-type antiferromagnetism, i.e. the spins of all 
nearest neighbors of a given lattice site belong to the 
other sublattice. This kind of antiferromagnetism is typ- 
ical in the strong coupling regime at and near half-filling 
because it allows for a maximum kinetic energy gain 
by virtual hopping processes, unlike a FM configuration 
which forbids these by Pauli's Principle. We assume the 
ground state has Neel structure, i.e. we set the two mag- 
netic sublattices to be saturated, {Sf) = S = —{Sf). 
Furthermore we do not take into account possible canted 
AFM configurations. 

We used the method proposed in Ref. to determine 
the FM/PS phase boundary.^^ This criterium for elec- 
tronic phase separation is based on the separation into 
AFM regions with one electron at each site and FM do- 
mains with an occupation Uc, a picture suggested by nu- 
merical results in Ref. On increasing n the AFM part 
grows until it occupies the whole system at half-filling. 
The total energy can be written as 



Utotin, w) = (1 - v)Uafm + vUfm 1 - 



1 



(23) 



where Utot is the total energy per site, Ufm is the FM 
energy per site and its argument is the particle density 
in the FM regions {n is the total particle density), Uafm 
is the AFM energy per site, and v = VFAi/Vtot is the 
FM volume fraction of the total system size. Minimizing 
the total energy with respect to v yields the following 
condition for the critical electron density ric at which 
electronic phase separation sets in (i.e. v = 1): 

Uafm = UpMiric) + (1 - nc)J7^M("-c) • (24) 

U'fm is ttie derivative of Ufm with respect to the par- 
ticle density. Note that we consider the electron density 
and not the hole density in the FKLM so that the corre- 
sponding formula in Ref. ^3 is modified accordingly. 

If one varies the chemical potential /i continuously the 
electron densities at which phase separation is present 
correspond to band occupations that cannot be stabilized 
as, e.g., demonstrated in the MC simulations in Ref. IT^ 
Given the jump in the electron density on varying fj, this 
kind of transition appears to be first-order. However, 
enforcing any value of n as in our case and having in 
mind the picture of AFM regions gradually taking over 
the whole system the transition from FM to AFM rather 
appears to be continuous. 
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FIG. 1: Total energy (in eV) for ferromagnetic, paramag- 
netic, and antiferromagnetic configuration for the square lat- 
tice. The AFM energy is evaluated at n = 1. Uc is computed 
according to the right hand side of Eq. 12411 . Parameters: 
S = 1/2, JS/t = 12. 



III. RESULTS AND DISCUSSION 

We computed phase diagrams for different spin quan- 
tum numbers S for a simple cubic lattice, a square lattice, 
and a ID chain at zero temperature. To simulate a classi- 
cal spin, i.e. a spin that can be oriented in any direction 
in space, we used a spin quantum number S = 10. In 
order to obtain a FM-PM-PS phase diagram we evalu- 
ated the total energy at T=0 for the saturated ferromag- 
netic, for the paramagnetic, and for the antiferromag- 
netic (n — 1) spin configuration of the core spins, respec- 
tively, as a function of the occupation number n and for 
several values of the Hund coupling J. To compare the 
results for different S we take the proper scaling cx JS of 
the interaction into account and normalize the localized 
spins. Thus in the following we consider localized spins 
S/S coupled to itinerant electrons by JS. 

Before starting the discussion of our results we should 
add two remarks. First, it is well known that ID systems 
exhibit some pecularities; for example, an integer spin 
nearest-neighbor Heisenberg chain has a gap in its exci- 
tation spectrum [Haldane gap). Moreover, non-local cor- 
relations are important, whereas our approach is based 
on a local self-energy. However, our approximate theory 
is applicable for any finite dimension and thus wc con- 
sider it worthwhile to present results for D = 1, too. 

Secondly, we have to address the issue of anisotropy. 
The Mermin- Wagner theoremSi forbids spontaneous 
symmetry breaking in _D < 3 for an isotropic Hamiltonian 
like ^ at finite temperatures. In order for our results to 
be relevant at small temperatures, too, we have to add 
an (infinitesimally) small anisotropy term, e.g. a single- 
ion anisotropy taking spin orbit coupling into account. 
Being orders of magnitude smaller than the leading en- 
ergy scale in our system, the Hund coupling J, it will 
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FIG. 2: PM (left) and FM (right) total energy (eV) of the 
square lattice for quantum 5 = 1/2 and classical localized 
spins. Parameters: JS/t — 24. 



not alter the phase boundaries visibly. There is another 
benefit of adding an anisotropy to On decreasing S 
the assumption of a Neel state for the antiferromagnetic 
phase becomes questionable due to quantum fluctuations. 
In ID this approximation even breaks down completely. 
These fluctuations are suppressed by anisotropy:-®^ 

Fig. Qshows the total energy of the FM and PM phases 
as a function of the band occupation n and of the AFM 
phase at n = 1. The result of the right hand side of Eq. 
(|^ is also plotted. The spin is 5* = 1/2 corresponding 
to a maximum of quantum fluctuations. The paramag- 
netic ground state at Ud = 0.45 emerges well before the 
criterium H24|) for PS is fulflUed at nc,2 = 0.72, indicating 
an instability of ferromagnetism against spin disorder. In 
our calculations wc did not find a second transition from 
PM to PS for 5 = 1/2 at n > nc,i. On the other hand for 
5 > 1 and the values of J we considered {JS/t > 6) we 
find that the critical value of n for the onset of phase sep- 
aration is always lower than the electron density where 
FM becomes unstable against PM, i.e. nc.2 < nc,i. 

To further analyse the FM-PM transition for S =1/2 
we show in Fig. [5] more results for the total energy and 
compare them to calculations based on classical localized 
spins. Whereas the PM energy is lower for the quantum 
spin over the whole range of electron densities the FM 
energies are practically the same for both spins. 

This can be related to the quasiparticle excitation spec- 
trum. As can be seen by Q the total energy of the 
FKLM is determined in complete formal analogy to the 
free electron case, i.e. by the (quasiparticle) density of 
states. Our findings suggest the following picture: In the 
FM phase there is a parallel alignment of the conduc- 
tion electrons with respect to the localized spins. Given 
a saturated FM spin background it is not possible for an 
itinerant electron to flip its spin by spin exchange. There 
is no signiflcant occupation in the |-band of the spectrum 
for any S. This is, at least for high Hund coupling, con- 
sistent with the excitation spectra we obtained (see FM 
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FIG. 3: Paramagnetic quasiparticle density of states of the 
square lattice. Parameters: JS/t=24, n=0.7. The inset shows 
the low energy band of the spectrum on a smaller scale. Ver- 
tical lines indicate the chemical potential. The dash-dotted 
line in the inset is the Bloch-shaped FM j-QDOS (5 = 1/2). 
The FM j-QDOS in this energy range practically vanishes 
(dotted line) indicating complete spin polarization. 

QDOS in the inset of Fig. |3J|. For the "f-electrons the lo- 
calized moments act as an effective field and their quasi- 
particle density of states (QDOS) is merely the Bloch 
density of states shifted by —^JS. Within this picture 
the scaling by JS is expected to transfer directly to the 
total energy which is indeed the case as can be seen by 
the practically identical curves in Fig. 121 (right). In other 
words, the scaling of the FM energy with JS expresses 
the fact that spin waves are frozen. 

The situation is different in the paramagnetic case. 
Now there are (energetically unfavored) states for t- and 
J,-electrons with an antiparallel orientation to the local- 
ized spins. Whereas the lower energy band in Fig. Ohas 
a "parallel character" of localized and itinerant magnetic 
moment the upper band corresponds to an antiparallel 
orientation. As there are more eigenstates with a paral- 
lel spin-spin alignment one expects larger spectral weight 
of the corresponding peaks in the excitation spectrum. 
The lower the magnitude of the localized spin the higher 
this difference becomes: the parallel case "outweighs" the 
antiparallel case most dominantly for S* = 1/2 ("triplet" 
vs. "singlet"). This can immediately be verified in the 
zero bandwidth limit. There is indeed a higher spec- 
tral weight for low spin quantum numbers as can be seen 
in Fig. 13 On the other hand increasing S results in 
an equal distribution of spectral weight for both bands 
in the classical limitiSS Thus, the paramagnetic state for 
lower magnitudes of the localized spins has lower total 
energy, as is observed in Fig. [21(left). 

Let us now proceed to the discussion of the magnetic 
phase diagrams we obtained in the spatial dimensions 
D =1,2, and 3. Fig. ^shows the phase diagram in _D = 
1. The phase boundaries indicate the transition from 
FM to the phase separated FM-AFM regime except for 



the case 5* = 1/2 which gives a FM-PM transition. For 
comparison we added results by other authors obtained 
by numerical methods. Our findings give the same overall 
picture as earlier works with an increasing FM region 
as the Hund coupling becomes stronger. However we 
find that the S — 1/2-FKLM has a significantly reduced 
ferromagnetic stability due to its "early" transition to 
a FM state compared to higher localized spin quantum 
numbers. A second important feature is that we do not 
see any major differences in our phase boundaries for 
S > 1/2. 

Our results compare well with the MC numerical phase 
diagram (classical spins) and with the DMRG results 
(5 = 3/2). There is no considerable variation of the 
phase boundary for S = 3/2 and S = oo either^ The 
deviations are largest for S = 1/2: here our theory ap- 
pears to underestimate the FM region. However, the 
authors of Ref. attribute the fact that their FM re- 
gion is reduced as compared to Ref. to the higher 
number of lattice sites they include in their computation. 
Hence it would be desirable to have more DMRG results 
with larger system sizes to see how the phase boundaries 
change. We did not find a second transition from FM to 
PS for S = 1/2, neither did we take spiral or incommen- 
surate correlations into account as was done in Ref. 
and Ref. However we can state that the reduction of 
the maximal region of FM for S = 1/2 with respect to 
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FIG. 4: Phase diagram in D = 1 for S = ^, | and classi- 
cal spin. Filled symbols refer to our results and indicate the 
FM/PS transition {S = 3/2, oo) and the FM/PM transition 
{S = 1/2). Open squares from Ref. ^3 mark the transition 
from FM to a spiral phase; in this work an island phase be- 
tween n « 0.2 and n « 0.8 and up to JS/t « 8, and PS 
for n > 0.8 and JS/t > 10 were also found. All other open 
symbols are taken from Ref. IT3 circles (5 = 1/2) and tri- 
angles {S = 3/2) indicate the boundary between FM and 
incommensurate correlations (IC), for n > 0.8 and above 
JS/t « 12 (5 = 1/2) and JS/t f« 6 (5 = 3/2) PS was 
found, too; diamonds {S = oo) mark transition from FM to 
PS {JS/t > 4) and from FM to IC at weaker coupling. Lines 
are guides to the eye. 
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FIG. 5: Phase diagram in Dj= 2 for S = |, | and classical 
spin. Open circles from Ref. {S = oo) indicate boundary 
between FM and PS above JS/t « 7 and between FM and 
a spiral phase at weaker coupling:. SB stands for Schwinger- 
boson. Open squares from Ref. [l3 {S = oo) mark transition 
from FM to PS above JS/t « 4 and from FM to IC at weaker 
coupling. Lines are guides to the eye. 



FIG. 6: Phase diagram in D = 3 for S = ^, | and classical 
spin. Open circles from Ref. indicate boundary between 
FM and PS above JS/t 7 and between FM and spiral phase 
at weaker coupling. Open squares from Ref. [ijj (semicircular 
density of states) mark transition from FM to PS. Lines are 
guides to the eye. 



higher spin quantum numbers is consistent with what we 
learn from the results obtained by other approaches. 

Fig. O and ini show the phase diagrams for 2D and 3D, 
respectively!^ They give essentially the same picture as 
in ID. Again we observe an unstability of the FM phase 
against PM for S = 1/2 only, reducing the region of 
FM stability as compared with higher S. There is no 
significant change of the phase boundary for 5* > 1/2 
apart from some enhancement of FM for classical core 
spins and larger Hund couplings in all dimensions. The 
S — oo results in 2D are in accordance with MC simu- 
lations in Ref. ^3 which yield FM for the full range of 
band filling at large J. We note a slight enlargement of 
the PS region in 3D for all S. However we do not want 
to emphasize the quantitative differences too much. As 
was already pointed out in Ref. Jdj the small differences 
in the energies of the different phases lead to uncertain- 
ties in the exact location of the phase boundaries. In our 
case we estimate these error bars to be about 10% with 
respect to the corresponding electron density n. For the 
same reason we are careful not to put too much signif- 
icance into the Z)-dependence of the crossing points of 
our 5 = 3/2 and classical S phase boundaries. How- 
ever it is interesting to note that there is a crossing in all 
dimensions. 

As before our findings agree well with those published 
by other authors who used different methods. There is 
one exception: compared to the other results the DMFT 
phase boundaries (however, for a, D = cx)-Bethe lattice) 
from Ref. ^ seem to overestimate FM considerably. 

We conclude that the magnetic properties of the 
FKLM at strong coupling and as far as the magnetic 
phases we investigated are concerned are rather insen- 



sitive to variations of the spatial dimension, at least at 
T=0. This falls in line with other results obtained using 
classical localized spins/i^ 



IV. SUMMARY AND OUTLOOK 

We have presented magnetic phase diagrams for the 
ferromagnetic Kondo lattice model in D — 1,2, and 3. 
Our method is based on an equation of motion decou- 
pling procedure fulfilling exact relations among Green 
functions and among their spectral moments. It does 
not require the assumption of classical spins. To deter- 
mine the phase boundaries we computed and compared 
the total energy of the different phases at zero tempera- 
ture. 

There are three main results. First the case S" = 1/2 
appears to be special exhibiting a reduced region of fer- 
romagnetism in the J-n-plane due to an instability of 
FM against spin disorder. Increasing the electron den- 
sity this transition always takes place before FM/AFM 
phase separation can occur. Secondly the phase bound- 
aries for S > 1/2 appear to be quite robust with respect 
to changes of the magnitude of the localized spin. This 
supports the widespread usage of classical localized spins 
in the treatment of the FKLM. Finally, these two features 
are recovered and quantitatively similar in the phase dia- 
grams in all dimensions we investigated, namely ID, 2D, 
and 3D. 

To our knowledge there are no numerical phase dia- 
gram results in 2D and 3D with quantum localized spins 
as this is numerically a quite demanding task. It would 
be interesting to explore if the same trends as in ID hold 
for S" = 1/2 using numerically exact methods like DMRG. 
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It should also be a worthwhile task to examine the influ- 
ence of the spin magnitude at higher temperatures up 
to the Curie temperature (in the vicinity of which colos- 
sal magnetoresistance occurs in the manganites). Phase 
boundaries may of course change to a certain extent with 
the crystal lattice structure, i.e. Bloch density of states. 
More changes can be expected when including a finite 
next-nearest neighbor hopping integral. Finally we fo- 
cussed on the phases thought to be relevant for the in- 
termediate to strong-coupling regime and left out other 
phases that come into play in the weak-coupling case. 
These issues are left for further investigation. 



Appendix A: Interpolation coefficients 

Exploiting spectral moment relations leads to the fol- 
lowing coefficients in the approximations (|17|I - H2U|) for the 
higher order local Green functions: 
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The mixed expectation values 
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can all be evaluated with the corresponding Green func- 
tions using the spectral theorem. 

Appendix B: Localized spin expectation values 



It holds for ferromagnetic saturation: 
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and for the paramagnetic phase: 
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